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SUMMARY 


An experimental investigation has been conducted which provides quantitative 
pressure and force data for five axisymmetric boattail nozzle configurations* These 
configurations simulate the variable-geometry feature of a single nozzle design oper- 
ating over a range of engine operating conditions. This investigation was performed 
at an angle of attack of 0° in the Langley 16-Foot Transonic Tunnel at Mach numbers 
from 0.60 to 1.30. The experimental data and theoretical predictions, provided by 
computational fluid- dynamics programs, were compared and found to be in generally 
good agreement. 


INTRODUCTION 

A supersonic cruise aircraft must be capable of operation over a wide altitude- 
velocity flight envelope, which includes such diversities as subsonic take-off and 
landing, subsonic cruise, climb and loiter, and supersonic cruise. In addition, 
maneuver requirements may be added for military aircraft. Obviously, current engine 
concepts cannot be optimized over such a broad range of requirements. One proposed 
solution for these diverse requirements is the variable-cycle engine (VCE) which 
incorporates a variable-engine geometry and combustion arrangement in order to oper- 
ate as a turbofan, an afterburning turbojet, a duct-burning turbofan, or an inter- 
mediate hybrid combination as required by any particular mission segment. Mission 
requirements, development concepts, and preliminary designs of the VCE may be found 
in references 1 to 8. 

The purpose of the current investigation was to provide quantitative pressure 
and force data for an axisymmetric variable-geometry nozzle which is applicable for 
use in the design of the engine for a supersonic cruise aircraft. This was accom- 
plished by testing five nozzles in the Langley 16-Foot Transonic Tunnel, which 
represent the variable-geometry feature of a single, convergent-divergent nozzle 
design over a range of engine operating conditions. Tests were conducted at nozzle 
pressure ratios from jet off to about 10 at free-stream Mach numbers from 0.60 to 
1.30. The experimental data were compared with theoretical predictions in order to 
help verify several existing computational flu id- dynamics programs. 


SYMBOLS 


e 

Aw 

A e 


increment of model cross-sectional area at metric break station 
67.31 cm, m 2 

o 

nozzle exit area, m 

maximum cross-sectional area of model, 182.415 cm 2 
nozzle geometric throat area, m 2 

2 

surface wetted area, m 

increment of axially projected boattail area, m 2 



^D,f 


- D, 0 


"F / i 


C 

P/ 


P 


boattail friction-drag coefficient (see section entitled "Data Reduction") 
boattail pressure-drag coefficient (see section entitled "Data Reduction") 
aerodynamic ideal thrust coefficient, F^/q^A^ 


boattail pressure coefficient. 



f / cyl 


e 

dm 


A, mom 


-bal 


F ± 


F j 


CO 

m 

%e 

P 

Pc 

Pe 

P t,j 


friction drag on cylindrical portion of model between stations 
67.3 1 cm and 137.16 cm, N 

nozzle total drag, N 

nozzle-exit diameter, cm 

model maximum diameter, 15.240 cm 

nozzle geometric throat diameter, cm 

thrust, N 

axial -momentum tare force, N 

axial force measured by balance, positive forward, N 
ideal isentropic gross thrust, N 
gross thrust, mV e + (p e - p 00 )A e , N 
axial length of boattail, cm 

axial length of nozzle convergent section, cm 

axial length of nozzle divergent section, cm 

free-stream Mach number 

nozzle measured mass-flow rate, kg/sec 

nozzle ideal mass-flow rate, kg/sec 

Rey no 1 ds numbe r 

local static pressure. Pa 

local static pressure in the metric break. Pa 
average static pressure in the nozzle-exit plane. Pa 
jet total pressure. Pa 
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(p t -i/Poo) design nozzle pressure ratio for ideally expanded 

^ es exhaust flow 

p local static pressure on the boattail, Pa 

P 

p^ free-stream static pressure. Pa 

free-stream dynamic pressure. Pa 

R gas constant for air, 287.3 J/kg-K 

r local radial distance from model center line, cm 

r model maximum radius, 7.62 cm 

m 

T. . jet total temperature, K 

T- 9 J 

V e average axial velocity in nozzle-exit plane, m/sec 

x axial coordinate with origin one model diameter upstream 

of the nozzle exit (positive downstream), cm 

P nozzle boattail angle, deg 

y ratio of specific heats, 1.3997 for air 

6 nozzle divergence angle (downstream of throat), deg 

9 nozzle convergence angle (upstream of throat), deg 

Abbreviations : 

A/B afterburning 

J/0 jet-off condition 

APPARATUS AND PROCEDURE 
Wind Tunnel 

The experimental investigation was conducted in the Langley 16-Foot Transonic 
Tunnel. This facility is a single-return, continuous -flow, exchange-air-cooled, 
atmospheric wind tunnel with an octagonal slotted-throat test section. It has a con- 
tinuously variable Mach number range from 0.20 to 1.30. A detailed description of 
this wind tunnel is given in reference 9. 

Models 

A set of five nozzle models was used to simulate the range of geometric 
settings of a variable-geometry axisymmetric nozzle applicable for supersonic cruise 
aircraft. A photograph of nozzle configuration 5 installed in the test section of 
the Langley 16-Foot Transonic Tunnel is shown in figure 1. Figure 2 depicts the 
general arrangement of the isolated nacelle model and support system. The five 
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nozzle configurations tested are shown in the photograph in figure 3. The simulated 
flight segment, simulated power setting, design nozzle pressure ratio, and design 
dimensions are presented in figure 4. 


Jet-Exhaust Simulation 

For jet-exhaust simulation, an external high-pressure air system provided a 
continuous flow of clean, dry air at a maximum pressure of 724 kPa and a controlled 
temperature of nominally 300 K ahead of the nozzle throat. As shown in figure 2, 
this high-pressure air was brought through the support sting and strut into a high- 
pressure plenum and was then introduced, through eight sonic nozzles, radially into 
the metric portion of the model to eliminate incoming axial momentum. Two flexible 
metal bellows provided an air seal between the metric and nonmetric portions of the 
nacelle. Finally, the air traveled through flow- smoothing screens into a stagnation 
chamber and expanded through the nozzle configuration being tested. 


Instrumentation 

Metric model forces and moments (excluding rolling moment) were measured by a 
five -component strain-gage balance. Jet total pressure and total temperature were 
measured by probes mounted in the metric tailpipe as shown in figure 2. Also, there 
was a row of internal static-pressure orifices from the nozzle entrance to the exit 
plane, and there was a row of external static-pressure orifices starting just forward 
of the nozzle boattail to the exit plane. The locations of the internal and external 
pressure orifices are given in table I. The pressure and temperature of the internal 
airflow were also measured in the high-pressure plenum located in the model forebody. 
(See fig. 2.) Static pressures in the gap at the metric break (station 67.31 cm) 
were obtained from 10 orifices distributed over the model cross-sectional area. 


Tests 

In accordance with the criteria of references 10 and 11, a boundary- layer tran- 
sition strip, 0.25 cm wide consisting of No. 100 silicon carbide grit sparsely dis- 
tributed in a lacquer film, was applied 2.54 cm downstream of the nose. Angle of 
attack was held at zero for all tests. Each nozzle configuration was tested at Mach 
numbers of 0.60, 0.80, 0.90, 0.94, 1.20, and nominally 1.30. The average Reynolds 
number per meter varied from 10.43 x 10 6 at a Mach number of 0.60 to 13.20 x 10 6 at a 
Mach number of 1.30. Nozzle pressure ratio was varied from jet off to approximately 
10 depending on the free-stream Mach number. 


Data Reduction 

All data for both the model and wind-tunnel facility were recorded simulta- 
neously on magnetic tape. The recorded data were used to compute standard force and 
pressure coefficients. All force coefficients in this report are referenced to the 
model maximum cross-sectional area. 

Nozzle boattail pressure drag was obtained during this test from pressures 
measured along the top of the boattail (see table I) by assigning an incremental, 
axially projected area to each orifice and using the stepwise summation equation 
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P B,k A (3,k 


( 1 ) 


°D,P 



A 

m 


where n is the number of orifices in the external row. To check these results, a 
number of continuous distributions of Cp ^ with (r/r m ) were manually integrated 
by planimeter and gave essentially the same drag values as equation ( 1 ) . Pressure 
forces on the small base at the nozzle exit are considered negligible and are not 
included in the boattail drag computation. 

Aerodynamic skin-friction drag coefficients were computed from the flat- 
plate formula for turbulent, compressible boundary layer given in reference 12. 


'D,f 


1 + 0.2M 


! )° 


0.472 

.467 


(log N ) 
10 Re 


2.58 


m 


( 2 ) 


At each test point, Cp f was calculated for Reynolds number N Re and wetted areas 
A^ corresponding to three different characteristic lengths on the model (fig. 2): 

(1) nose to exit, (2) nose to the upstream end of the boattail at station 137.16 cm, 
and (3) nose to metric-break station 67.31 cm. The difference between calculations 
(1) and (2) provided the estimated skin-friction drag of the boattails; and the dif- 
ference between calculations (2) and (3) was used to obtain the friction drag force 
on the cylindrical part of the metric afterbody D f As shown later, cyl 

was used to correct the balance data. 

Actual mass -flow rates used to compute discharge coefficients m/m. and ideal 
thrust for the test nozzles were measured by the system of eight radial nozzles 
supplying air to the metric afterbody. This required a pretest calibration of the 
internal-flow system which will be discussed subsequently. The ideal mass-flow rates 
for the test configurations were computed from stagnation pressure and temperatures 
measured in the tailpipe by using the choked-flow equation 



Ideal thrust, defined as the product of measured mass flow and the ideal isen- 
tropic velocity, was obtained from the relation 



(4) 


5 


Since the model strain-gage balance (fig. 2) measures the sum of pressure 
and viscous forces on the entire metric afterbody (model portions aft of station 
67.31 cm), a number of corrections must be applied to the raw-balance data in order 
to isolate the sum of external and internal forces on the nozzle alone. Nozzle gross 
thrust-minus -drag performance was computed from the relationship 


F . 
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10 


= F. 


bal 


(P 


k=1 


c,k 


- p )A 


c,k 


+ D 


f ,cyl 


- F 


A, mom 


(5) 


In equation (5) the term represents the raw-balance output corrected for 

interactions and model weight tares. The pressure-area term corrects for the force 
on the front face of the metric afterbody caused by differences between the cavity 
pressure p c in the metric break and free-stream pressure p^* Even though experi- 
ence with this model has shown cavity pressures to be extremely uniform, 10 orifices 
were used in these tests to determine this balance-correction term. As indicated 
previously, the term cyl is the friction drag on the cylindrical section of the 

metric afterbody, which must be restored to the balance reading since this force is 
not associated with the nozzle. The term F^ mom corrects for axial momentum and 
bellows tare forces caused by the high-pressure air flowing from the nonmetric high- 
pressure plenum into the metric tailpipe through the eight small injection nozzles. 
(See fig. 2.) Although the air is injected radially and the flexible seals (metal 
bellows) are placed in tandem in an effort to eliminate such tares, small forces do 
arise in practice and must be taken into account. In the present investigation these 
forces were generally less than 2 percent of ideal thrust. 

Axial-momentum tares were evaluated by statically (M^ = 0) testing several 
"standard” exhaust nozzles over the range of supply pressures, mass-flow rates, and 
throat areas required by the five test nozzles. The "standard" nozzles were inter- 
nally convergent, with a circular-arc longitudinal profile tangent to the axial 
direction at the exit, and had a profile arc radius equal to twice the exit 
diameter. (See ref. 13.) The correct values of Fj/F^ and m/nu as functions of 
nozzle pressure ratio have been well established for these "standard" nozzles from 
previous tests and were used to predict gross thrust for comparison with the balance 
measurements. Any force difference, or momentum tare, was computed from 
equation (5), which for static conditions became simply 


F = F - F 

A, mom bal j 


( 6 ) 


Resulting tare forces were correlated as a function of pressure measured in the high- 
pressure plenum in the model nose, thus providing a calibration from which tares were 
computed for the test-nozzle data. 

Internal-flow measurements taken during the aforementioned tare calibrations 
were also used to compute discharge coefficients for the system of eight radial 
nozzles supplying air to the tailpipe. The actual mass-flow rates were those 
measured by the "standard" nozzles installed at the rear of the model. Since the 
eight air-supply nozzles were choked at all test conditions, the ideal mass-flow 
rates were computed from equation (3), by using the total throat area, and from the 
high-pressure plenum measurements as stagnation conditions. Resulting discharge 
coefficients were then plotted against plenum pressure, thus providing a calibration 
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of the air-supply system as a sonic flow meter. This calibration, and measurements 
taken in the high-pressure plenum, were later used to measure airflow when investi- 
gating the five test nozzle configurations. 

Following the aforementioned calibrations, the five test nozzles were investi- 
gated at static (M^ = 0) conditions, and equation (6) was used to compute gross 
thrust. Wind-on gross thrust was also computed by using the wind-tunnel data and 
combining results of equations (1), (2), (4), and (5) in the equation 


F . 
3 

F, 

r 



c + c 
D,P D,f 


F,i 


(7) 


where 


~F,i 


is the aerodynamic ideal thrust coefficient, 


PRESENTATION OF EXPERIMENTAL AND ANALYTICAL RESULTS 

The data from the experimental tests are presented in figures 5 to 18. Measured 
pressure distributions are shown in figures 5 to 9. Various measured force param- 
eters are presented in figures 10 to 18. Results from computational fluid- dynamics 
programs are compared with selected experimental data in figures 19 to 25. A data 
list of the figures is presented as follows: 


Figure 


External and internal surface-pressure distributions 

of nozzle configuration 1 5 

External and internal surface-pressure distributions 

of nozzle configuration 2 6 

External and internal surface-pressure distributions 

of nozzle configuration 3 7 

External and internal surface-pressure distributions 

of nozzle configuration 4 8 

External and internal surface-pressure distributions 

of nozzle configuration 5 9 

Variation of jet-off, boattail pressure-drag coefficient with 

free-stream Mach number 10 

Influence of nozzle pressure ratio on boattail pressure-drag 

coefficient 11 

Calculated skin-friction drag of the various boattails 12 

Nozzle internal performance at static and wind-on conditions 13 

Static discharge coefficients of the nozzles 14 

Measured thrust-minus -drag performance of all nozzle-boattail 

combinations 15 

Aerodynamic ideal thrust coefficient as a function of jet 

total-pressure ratio 16 

Representative schedule of nozzle pressure ratio with Mach number 
for a low-bypass -ratio turbofan engine at maximum rotational 

speed •••• 17 

Summary of nozzle-boattail overall performance for the 

pressure-ratio schedule of figure 17 •••••••••••••••• • • • 18 

Comparison of experimentally and analytically determined 

external pressure distributions for nozzle configuration 1 19 
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Figure 


Comparison of experimentally and analytically determined external 

pressure distributions for nozzle configuration 2 20 

Comparison of experimentally and analytically determined 

external pressure distributions for nozzle configuration 3 21 

Comparison of experimentally and analytically determined external 

pressure distributions for nozzle configuration 4 22 

Comparison of experimentally and analytically determined external 

pressure distributions for nozzle configuration 5 23 

Comparison of experimentally and analytically determined internal 

pressure distributions 24 

Comparison of experimentally and analytically determined variation 

of static thrust coefficient with nozzle pressure ratio 25 


DISCUSSION OF EXPERIMENTAL RESULTS 
Pressure Measurements 

The external and internal surface-pressure distributions are given for each of 
the five nozzle configurations in figures 5 to 9. The nozzle longitudinal cross- 
section contour has been included in each figure to show the starting point (station 
137.16 cm) and angle of the boattail for the external pressure distributions and to 
show the throat location and divergence angle for the internal pressure distribu- 
tions • The symbol J/O in the keys of figures 5 to 9 denotes the jet-off condition 
and pertains to external pressures only. 

External pressur e distributions . - Jet operation generally had a favorable effect 
on the external pressure coefficients for the nozzles with low expansion ratios 
(A e /A t = 1.25 to 1.50) and large boattail angles from 8° to 15° (configurations 1 , 

3, and 4). Results for these configurations at M^ = 0.90 and 0.94 (for example, 
see fig. 5(c)) also indicate a standing shock wave located between x/d m = 0.2 
and 0.3, depending on configuration, M^, and p t ./p^. Increasing p t j/p^ 
tends to push the shock wave upstream, and external-flow separation occurs downstream 
of the shock wave for these configurations. At supersonic speeds (see fig. 5(e)), a 
trailing-edge shock occurs on the nozzle near the exit. As p t j/p^, is increased, 
the trailing-edge shock is pushed upstream and large jet interference effects are 
observed behind the shock. Jet operation had little or no effect on the forward 
portion of the boattail pressure distributions at supersonic speeds, because the 
downstream disturbances can only feed forward through the subsonic boundary layer. 

Separation of the external flow is shown to be primarily dependent upon boattail 
angle, because the data for configurations 2 and 5 with (3 < 4° showed no evidence 
of shock formation or flow separation at any test condition. With increasing 
P t j/Pa>' boattail surface pressures decrease generally, indicating unfavorable 

jet effects. This is attributed to strong turbulent mixing and entrainment in the 
vicinity of the exit resulting from overexpansion and separation of the jet flow 
internally. Jet-interference effects extended well forward on the boattails of con- 
figurations 2 and 5 at subsonic speeds, but they were confined to the last pressure 
orifice just ahead of the exit at M^ = 1.20 and 1.28. Thus, the expansion occurring 
along the rear of the boattail at subsonic speeds gradually concentrates into an 
expansion fan centered on the nozzle-exit lip at supersonic Mach numbers. 

Internal pressure d istri butions .- The ratio of static pressure to total pres- 
sure measured along the inside surface of the nozzles is shown at the bottom of 
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figures 5 to 9. Corresponding values for uniform / one-dimensional flow were computed 
from the the gas tables of reference 14 for the internal local-to-sonic area ratios 
of configuration 1. These results are compared with the experimental data in fig- 
ure 5(a) and show that the axisymmetric flow overexpands noticeably at the throat 
(x/d m = 0,2)* This is generally characteristic of nozzles with relatively steep 
convergent sections and "sharp" throats (small profile radius). Previous analytical 
studies, such as reference 15, have shown that the flow tends to "overturn" when 
passing around the throat profile. Therefore, the local flow angles relative to the 
nozzle axis exceed the wall divergence angle in the region behind the throat near the 
wall. Compressive turning is then required to match flow and wall angles, thus 
resulting in the positive pressure gradient behind the throat evident in much of the 
experimental data. Compression waves from this region coalesce in the vicinity of 
the nozzle axis and tend to form an internal shock surface that either intersects the 
wall further downstream or passes through the exit opening, depending on nozzle 
geometry. Evidence of the disturbance reaching the wall can be seen in figure 5 at 
x/d^ = 0.7 and in figure 7 at x/d m “ 0*85 (configurations 1 and 3, respec- 
tively). Notice that the small pressure rise at each of these locations appears only 
when the internal flow is locally attached? it is not, therefore, related to separa- 
tion from overexpansion. 

As the divergence angle increases, thereby increasing the flow area more rapidly 

behind the throat, the walls are better able to accommodate the "overturning" of the 

flow. The wall pressure gradient leaving the throat decreases correspondingly, and 
the tendency of this region to form internal shocks disappears. (For example, 
compare figs. 5 and 6.) 

The internal pressure distributions show shock-induced flow in the divergent 
section of the nozzle for nozzle pressure ratios significantly less than the design 
value. Increasing the divergence angle 6 of the nozzle not only increases the 

axially projected area of the divergent section but also moves the separation loca- 
tion forward at a given and nozzle pressure ratio. The latter was partic- 

ularly noticeable as 5 varied from 4.8° to 13.2° (configurations 4, 5, and 2), 
and this effect may be seen by comparing corresponding parts of figures 8, 9, and 6, 
in that order. Thus, increased divergence and extensive separation combine to expose 
increasing amounts of internal surface area to local pressures surrounding the exit. 
This, in turn, causes the external flow to have stronger effects on the pressure 
forces generated by the divergent section. 

The largest effects of on internal pressures occurred with configura- 

tions 2 and 5, which were the most extensively separated nozzles. Figures 6 and 9 
show that increasing at fixed nozzle pressure ratio reduced the pressure level 

in the separated region and tended to move the separation location rearward. Also, 
the most noticeable pressure reduction in the separated region occurred as 
increased from 0.94 to 1.20 and the external flow began to undergo supersonic expan- 
sion at the exit lip. Consequently, increasing can be expected to increase the 

overexpansion and drag associated with the divergent section at low nozzle pressure 
ratios. 

Obviously, as nozzle pressure ratio is increased by increasing the jet stag- 
nation pressure p t j, the separation moves to the exit and disappears. Once the 
exit flow is fully supersonic, the internal pressure distribution becomes independent 
of external conditions. For the present configurations, the minimum nozzle pressure 
ratio required to eliminate separation ranged from approximately 3 to 7 for corre- 
sponding design values of 4.25 to 21.23 and for expansion ratios of 1.25 to 3.00, 
respectively • 
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External Performance 


The variation of jet-off, boattail pressure-drag coefficient with free-stream 
Mach number is shown in figure 10. Jet-off boattail pressure drag generally 
increases with increasing nozzle boattail angle. As the Mach number approached 
unity, the steepest drag rise occurred for nozzles with p = 8° to 15°, starting 
at » 0.80. As might be expected, the shallow boattails of configurations 2 

and 5 exhibit a much smaller drag rise, which started at M^ » 0.90 or slightly 
greater. In the supersonic-speed range, boattail pressure-drag coefficient is almost 
directly proportional to nozzle boattail angle. 

Figure 1 1 presents the influence of nozzle pressure ratio on boattail pressure- 
drag coefficient for all configurations and test Mach numbers. In general, these 
curves are shaped by two opposing effects: the growth in jet diameter, or pluming 

effect, which pressurizes the boattail surface and thus reduces drag? and jet 
entrainment, which tends to accelerate the flow on the boattail, thereby increasing 
the drag. Therefore, where the slopes of the curves are negative, the pluming effect 
predominates, and where positive, the influence of jet entrainment predominates. On 
any given curve, a drag level less than the initial jet-off value (p t j/Pco ** ^ 
indicates a favorable jet effect. Furthermore, all subsequent comments regarding jet 
effects on boattail drag refer to the portion of the drag curves between p t j/Poo = 3 
and the maximum test value, since this is the range of practical interest for the 
Mach numbers of this test. 

Jet operation had a favorable effect on the boattail drag of nozzle configura- 
tions 1, 3, and 4 at all test Mach numbers. At subsonic conditions the jet effects 
on drag are caused by modest increases in pressure which extend over all or most of 
the boattail length at a given nozzle pressure ratio. Among configurations 1, 3, 
and 4, the jet effects at M^ = 0.60 and 0.80 are seen to increase with increas- 
ing boattail angle (projected area) such that the drag levels are nearly equal for 
P = 8° to 15°. As M^ becomes supersonic, the influence of the jet on boattail 
pressures becomes limited to the rear half or less of the boattail length. However, 
as the nozzle pressure ratio increases, the pluming jet pushes the trailing shock 
forward, thus substantially increasing the pressure and the area affected behind the 
shock. Consequently, the effect of the jet on the drag coefficients of configura- 
tions 1, 3, and 4 is still strongly favorable at supersonic speeds as reflected in 
figure 11(c). 

The shallow boattails of configurations 2 and 5 exhibit much smaller variations 
in pressure-drag coefficients with nozzle pressure ratio than the other configura- 
tions. At subsonic speeds, jet effects on the drag of configurations 2 and 5 are 
unfavorable. The increase in boattail drag with jet operation is attributed to the 
dominance of jet entrainment when the nozzles are operating highly overexpanded 
(P t j/Poo <<: Design value) and, therefore, produce little, if any, beneficial pluming 
effect. At supersonic speeds the boattail pressure data for configurations 2 and 5 
showed no shock formation or separation, and jet effects were confined to a very 
small area just ahead of the exit. Figure 11(c) confirms that jet effects on drag 
at M^ = 1.20 and 1.28 are negligible. 


Internal Performance 

Nozzle internal thrust performance F j/ F i was determined by equation (7). 
Thrust was obtained by combining the boattail pressure drag of figure 11 and boat- 
tail skin-friction drag shown in figure 12 with the thrust -minus -drag measurements 
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from the model balance. The resulting variations of thrust ratio with nozzle pres- 
sure ratio are indicated in figure 13, and the discharge coefficients are shown in 
figure 14. The internal performance of the test configurations may be summarized as 
follows: 

(1) Nozzle configurations 1 , 3, 4, and 5 achieved an average peak thrust per- 
formance of F j/ F i » 0.99 at values of nozzle pressure ratio at or near design value 
(optimum expansion). 

(2) Experimental nozzle pressure ratios for nozzle configuration 2 were substan- 
tially lower than the design value; therefore, this configuration suffered overexpan- 
sion losses throughout the range of nozzle pressure ratios shown. However, it is 
noted that Fj/F^ « 0.95 at 50 percent of design nozzle pressure ratio, which is 
characteristic of all configurations when free of external-flow effects. (See con- 
figurations 1, 3, 4, and 5 at M = 0.) Therefore, the peak performance of configu- 
ration 2 would also be expected to reach F j/ F i « 0.99 at the design point. 

(3) When the internal flow was considerably overexpanded and separated, the 
presence of external flow tended to reduce the thrust performance of the nozzles 
relative to the wind-off case (M^ = 0), particularly at supersonic speeds. At a 
given nozzle pressure ratio, this additional overexpansion loss increased with 
increasing divergence angle and free-stream Mach number as indicated by the internal 
pressure data discussed previously. 

(4) At nozzle pressure ratios high enough to eliminate internal separation, the 
thrust ratios for all Mach numbers converged to one curve. This indicated fully 
supersonic exit flow with internal performance independent of external conditions. 
Reference 16 shows similar effects of internal-flow separation for another group of 
axisymmetric, convergent-divergent nozzles. 

(5) Figure 14 shows that the discharge coefficients of all nozzles are essenti- 
ally independent of nozzle pressure ratio (and throat Reynolds number) when operating 
choked. Values of m/m^ ranged from 0.964 (configuration 1) to 0.976 (configura- 
tion 5), the average level being rather typical of nozzles with relatively "sharp" 
throats. The improvement in discharge coefficient between configurations 1 to 5 is 
attributed to more uniform flow entering the throat as the convergence angle 
decreased and also to the boundary layer thinning as the pressure gradient behind the 
throat decreased with increasing divergence angle. 


Overall Performance 

The ratio of nozzle thrust minus drag to ideal thrust plotted against nozzle 
pressure ratio for all configurations and test Mach numbers is shown in figure 15. 

The shape of these curves is determined primarily by thrust characteristics, and the 
maximum overall force ratio approaches the internal performance Fj/F^ as F^ 
becomes large relative to drag at high pressure ratios. At subsonic speeds, the 
relative levels of the data reflect the nearly equal drag characteristics of configu- 
rations 1, 3, and 4 (fig. 11(a)) as well as the large overexpansion losses of config- 
urations 2 and 5 (fig. 13), both cited previously. As the free stream becomes super- 
sonic, higher drag losses and the inverse variation of overall performance with 
boattail angle become increasingly evident. 


Aerodynamic ideal-thrust coefficient of all nozzle configurations tested is 
plotted in figure 16 as a function of nozzle pressure ratio for the various test Mach 
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numbers. Since configurations 1 and 2 have identical throat areas and their static 
discharge coefficients are nearly the same (see fig. 14), their ideal thrust coeffi- 
cients, within the graphical accuracy of figure 16, appear as one plot. (See 
fig. 16(a).) These graphs provide the means of converting aerodynamic coeffi- 
cients based on q co A m to propulsion-force ratios based on F^. (That is, 

F /^ m C F,i = F /*V> 

Figure 17 presents a representative schedule of nozzle pressure ratio with Mach 
number for a low bypass-ratio turbofan engine operating at maximum rotational speed 
(climb and acceleration conditions). This schedule was obtained from unpublished 
industry sources used in performance studies of advanced supersonic aircraft. Fig- 
ure 18 shows the overall performance (fig. 16) for the nozzle configurations 
appropriate to the test speed range at typical operating pressure ratios and Mach 
numbers of the schedule (fig. 17). For convenient reference, the corresponding 
value of nozzle pressure ratio from the schedule is repeated with each Mach number 
at the bottom of figure 18. Also indicated are the optimum expansion ratios obtained 
from reference 14 for the schedule pressure ratios at each end of the subsonic and 
supersonic Mach number segments. The solid symbols denote internal performance for 
configuration 5 and show the rate at which the thrust optimizes for A e /A t = 2.0 
as schedule conditions increase toward the design point (p t j/p^ = 10.64) which is 
reached near M^ = 1.6. (See fig. 17.) By using figure 18, the overall performance 
can be summarized for a variable -geometry nozzle operating at optimum expansion 
(Fj/Fi = 0.99) along the assumed flight schedule of figure 17. 

At = 0.60 to 0.80, configuration 1 is operating at optimum internal expan- 

sion (F j/F^ = 0.99), which can be seen in figure 13 at schedule pressure ratios. 

If the exit areas of configurations 3, 4, and 5 were adjusted to provide the same 
internal conditions, it can be shown that the resulting boattail angles of all 
configurations then range from approximately 8° to 15°. It was shown previously 
(fig. 11(a)) that boattails in this angle range have nearly equal drag values at 
M ffi = 0.60 and 0.80. Therefore, all optimized configurations would be expected to 
have the same performance loss due to boattail drag as that shown by the data for 
configuration 1. As a result, the overall force ratio for optimum expansion at 
all power settings is approximately 0.98 as shown for configuration 1 between 
M =0.60 and 0.80. 

oo 

At Mach numbers between 0.80 and 0.94, test configurations 1, 3, and 4 are all 
operating at or near optimum expansion. (See fig. 13.) Therefore, the experimental 
data curves show overall force ratios for this condition directly, as well as the 
performance loss from boattail drag alone as each configuration enters the drag 
rise. It can be shown from the nozzle geometry and drag data that optimizing the 
exit area of configuration 5 for this Mach number segment results in a boattail 
angle and overall performance very similar to that of configuration 4. The flagged 
symbol at M^ = 0.90 denotes a representative subsonic cruise point at reduced power 
(i.e., not associated with the schedule). The cruise pressure ratio assumed was 
Pt j/Poo “ 4/ is an average of values generally associated with advanced 

fighters and transports. At = 0.90, overall force ratios are shown to range 

approximately from 0.95 at reduced-power cruise to 0.975 for afterburning conditions. 

At supersonic speeds, the overall performance for optimum expansion was esti- 
mated from the test data. For each Mach number and throat area, boattail angles 
were determined from the nozzle geometry for the optimum- expans ion ratios of 1.58 
and 1.69 required by the schedule. From cross plots of drag against boattail angle, 
the values of boattail drag for optimum- expans ion settings were selected. These 
values of boattail drag were converted to D n / F j_ by dividing them by the appropriate 
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values of C F ^ (fig. 16) and then subtracting the result from F./F^ = 0.99. The 
results, labeled "Optimum expansion" in figure 18, are identified by line code with 
the corresponding throat area of the test configurations with fixed-expansion ratio. 

For A t /A m = 0.42, the experimentally measured overall performance for 
A e /A t = 2.00 is seen to equal that estimated for optimum- expans ion ratios of the 
schedule. This indicates that the rates of change of thrust and drag with expansion 
ratio have equalized at these speeds. Consequently, if the exit area is increased 
somewhat beyond the optimum value, the resulting internal overexpansion loss is off- 
set by an equal reduction in boattail drag. Moreover, the calculations verify that 
this is true for all three throat areas. At M^ = 1.3 it can be seen that the over- 
all force ratio for optimally expanded settings from maximum dry power to maximum 
afterburning should range from 0.96 to 0.98, respectively. Based on the experimental 
data at = 1.3, configuration 5 should provide an overall force ratio of approxi- 

mately 0.985 when the schedule and design pressure ratios match in the vicinity of 
= 1.6. (See fig. 17.) This is verified in reference 16 by performance data for 
similar convergent -diver gent nozzles tested at higher supersonic Mach numbers. In 
conclusion, figure 18 shows the general importance of maintaining optimum expansion 
in flight and the high performance that results over a wide range of power settings. 


ANALYTICAL INVESTIGATION 

The following computer programs were used to provide theoretical comparisons 
with the experimental data: 

DONBOL (Douglas-Neumann/Boundary Layer) : This program, described in refer- 

ence 17, was used to calculate external surface-pressure distributions over the axi- 
symmetric boattail nozzles in the subsonic-speed regime (M^ = 0.60 and 0.80). The 
DONBOL program couples a Neumann solution for inviscid external flow, a modified 
Reshotko-Tucker integral boundary- layer technique, the control- volume method of Presz 
for calculating flow in the separated region, and an inviscid one-dimensional solu- 
tion for the jet-exhaust flow. The viscous and inviscid flows are solved iteratively 
until convergence is obtained. This method is limited to subsonic free-stream Mach 
numbers below that for which the flow over the body of revolution becomes sonic; 
also, this method does not include the effects of jet entrainment. Specifying the 
separation location is not necessary; however, the start and end of the region 
searched for separation must be specified. 

MOC (Method of Characteristics): This program was used to calculate the after- 

body surface-pressure distribution in the supersonic-speed regime (M^ =1.2 
and 1.28). The MOC program performs an axisymmetric method-of-characteristics 
solution for supersonic flow over boattailed afterbodies and is restricted to axi- 
symmetric, inviscid, irrotational flow. Since a description of this program, which 
was developed by Lawrence E. Putnam, Bobby L. Berrier, and Richard G. Wilmoth of the 
Langley Research Center, has not been published, reference 18 is given as a descrip- 
tion of the theoretical development of the axisymmetric method of characteristics. 

NAP (Nozzle Analysis Program): This program, described in reference 19, was 

used to calculate the nozzle internal pressure distribution and gross internal 
thrust. The NAP is capable of calculating time -dependent, perfect gas, inviscid, 
steady, and unsteady flow in two-dimensional and axisymmetric nozzles by using the 
nonconservative form of the two-dimensional equations of motion for a geometric 
computational grid which spans the flow area. Interior mesh points are calculated by 
using a MacCormack finite-difference method (a two-step predictor-corrector method 



described in ref, 20); a characteristics method is used to compute the boundary mesh 
points; and shock waves are modeled by using a shock-smearing technique which incor- 
porates an explicit artificial viscosity term with numerical smoothing described in 
reference 21. Gross thrust is obtained by integrating the momentum flux and dif- 
ference between exit and ambient pressures over the exit area. The NAP becomes 
increasingly inaccurate below design nozzle pressure ratio since limitations in the 
code prevent the flow velocity from returning to subsonic levels once it has become 
supersonic. This investigation is an example of NAP being applied to the analysis of 
an axisymmetric converging-diverging nozzle. Reference 22 gives an example of NAP 
being applied to the analysis of a nonaxi symmetric converging- diverging wedge nozzle. 


DISCUSSION OF ANALYTICAL RESULTS 
External Pressure Distribution 

A comparison of experimentally and analytically determined external pressure 
distributions is given in figures 19 to 23 which present data for nozzle configu- 
rations 1 to 5/ respectively. For part (a) of each figure which presents free-stream 
Mach numbers of 0.60 and 0.80, the theoretical prediction was computed by the DONBOL 
program of reference 17. Reasonable agreement with the experimental data trends was 
obtained for nozzle configurations with the larger boattail angles ( |3 > 8°). The 
DONBOL program generally underpredicted the strength of the nozzle shoulder expansion 
and the compression near the nozzle exit. Poor agreement between theory and 
experiment was obtained for the two nozzles with very small boattail angles 
(p < 4°). This discrepancy between the predictions of the DONBOL code and experiment 
for the small boattai 1-angle nozzles (configurations 2 and 5) is associated primarily 
with the analytical modeling of the jet-exhaust plume. The simple one-dimensional- 
solution technique for the jet- exhaust plume used in the DONBOL code is not adequate 
for the highly overexpanded nozzle cases of the present investigation. Also, the 
neglect of jet-entrainment effects by the DONBOL code is a factor in the observed 
discrepancies especially near the nozzle exit. Part (b) of figures 19 to 23 presents 
a conparison of theory and experiment for free-stream Mach numbers of 1.20 and 1.28 
for configurations 1 to 3 (1.20 and 1.27 for configurations 4 and 5). 

Good agreement was obtained for configurations with small boattail angles (con- 
figurations 2 and 5) except near the nozzle exit. For configurations with large 
boattail angles (configurations 1, 3, and 4), good agreement of theory and experiment 
was obtained until a shock wave and/or separated region was encountered on the nozzle 
external surface. Large discrepancies were noted in these regions. It is noted that 
the MOC computer code utilized is an isentropic, inviscid method-of-characteristics 
program. Therefore, although this program can reasonably predict the surface- 
pressure distribution of configurations with small boattail angles, it cannot predict 
surface pressures behind shock waves or in separated-f low regions associated with 
nozzle configurations having large boattail angles. 


Internal Pressure Distribution 

A comparison of experimentally and analytically determined internal pressure 
distribution for each nozzle configuration is given in figure 24. The experimental 
data shown were all taken at a free-stream Mach number of 0.60. However, internal 
pressure distributions are independent of free-stream Mach number at and above design 
nozzle pressure ratio. Below design nozzle pressure ratio, the experimental data 
may show effects of flow separation which cannot be predicted by the inviscid NAP of 
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reference 19. Good agreement is shown in figure 24 between experiment and theory, 
except in the converging section just upstream of the throat of each nozzle. In 
this section of the nozzle, the analytical predictions are consistently low. For 
configuration 2 which has a design nozzle pressure ratio of 21.23, the agreement 
between experiment and theory was remarkably better than expected, since the experi- 
mental nozzle pressure ratio of 6.97 was significantly lower than the design value. 
However, the nozzle was flowing full up to approximately x/d m = 0.95, thereby making 
the effects of separation minimal. 


Static Thrust Coefficient 

Figure 25 shows a comparison of experimentally and analytically determined vari- 
ation of static thrust coefficient with nozzle pressure ratio for each nozzle con- 
figuration. The solid line shown as theory represents the internal gross thrust 
computed by the NAP of reference 19. The dashed line depicted as ideal thrust was 
calculated by the equation for optimum thrust coefficient 





( 8 ) 


As expected, the experimental and theoretical values were less than the ideal 
value for all configurations. Excellent agreement between theory and experiment was 
obtained for nozzle configurations with low expansion ratios (configurations 1, 3, 
and 4). However, it is noted that for configurations 2 and 5 (high expansion 
ratios), the experimental values significantly exceeded the theoretical values at the 
lower nozzle pressure ratios tested. This may appear at first as a data anomaly when 
considering that the NAP does not consider viscous losses. The anomaly is best 
understood by considering that the maximum nozzle pressure ratio tested for these 
configurations is far below the design values of 21.23 and 10.69 for configurations 2 
and 5, respectively. 


Examination of the experimental internal pressure distribution for these two 
nozzle configurations (figs. 6 and 9 for configurations 2 and 5, respectively) 
readily shows the large amount of separation present. Large shock-induced separated- 
flow regions are indicated for these configurations at nozzle pressure ratios below 
7.0 for configuration 2 and below 4.0 for configuration 5. Since the NAP is limited 
to inviscid flow, and flow velocities cannot return to subsonic levels once super- 
sonic, NAP obviously cannot provide accurate results when large separated regions 
exist. However, the trends shown in figure 25 do indicate that NAP is capable of 
providing accurate results for these nozzle configurations at higher nozzle pressure 
ratios. Also, reasonable predictions far below design nozzle pressure ratio are pos- 
sible. This justifies using the code over a wide range of nozzle pressure ratios 
during preliminary design and development of new nozzles- 


CONCLUSIONS 

An experimental investigation has been conducted which provides quantitative 
pressure and force data for five axisymmetric convergent- diver gent nozzles with dif- 
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ferent internal and external geometry. These nozzles simulate the variable-geometry 
nozzle operating over a range of engine operating conditions. This investigation was 
conducted at an angle of attack of 0° in the Langley 16-Foot Transonic Tunnel at 
free-stream Mach numbers M^ from 0.60 to 1.30. Nozzle pressure ratio was varied 
from jet off to 10, depending on configuration and Mach number. The experimental 
data were also compared with computational theoretical predictions for internal pres- 
sure distribution, gross internal thrust, and external pressure distribution at sub- 
sonic and supersonic free-stream Mach numbers. 

The results of this experimental and analytical investigation indicated the 
following conclusions: 

1. Jet-off, boattail pressure-drag coefficient was generally proportional to the 
nozzle boattail angle. 

2. The amount in external-flow separation on the nozzle generally increased with 
increasing boattail angle, free-stream Mach number, and nozzle pressure ratio. 

3. Configurations which have very low boattail angles show only minor variation 
of pressure drag with nozzle pressure ratio, especially at supersonic speeds. 

4. Separation of internal flow was directly proportional to divergence angle and 
inversely proportional to nozzle pressure ratio. 

5. The ratio of thrust minus drag to ideal thrust for optimum expansion at 

Moo = was 0-95 for subsonic cruise and 0.975 for afterburning throat settings. 

At = 1.3, the overall force ratio for optimum expansion varied from 0.96 to 0.98 

for throat areas from maximum dry power to maximum afterburning, respectively. 

6. At subsonic speeds, theoretical external pressure distributions predicted by 
the DONBOL (Douglas-Neumann/Boundary Layer) program were in reasonable agreement with 
the experimental data for configurations which had boattail angles greater than 8° 
but were in poor agreement with experimental data for highly overexpanded nozzle 
configurations with boattail angles less than 4°. 

7. At supersonic speeds, theoretical external pressure distributions predicted 
by the MOC (Method of Characteristics) program were in good agreement with the 
experimental data until a shock wave and/or separated region was encountered on the 
nozzle external surface. 

8. Theoretical internal pressure distributions predicted by the NAP (Nozzle 
Analysis Program) were in good agreement with the experimental data except in the 
converging section (upstream of throat) of each nozzle. In this section of the 
nozzle, the analytical predictions are consistently low. 

9. The theoretical static thrust coefficient, predicted by the NAP, was in 
excellent agreement with the experimental data except at nozzle pressure ratios which 
produce large internal-flow-separation regions. The nozzle pressure ratios which 
produce internal-flow separation were found to be well below the design nozzle pres- 
sure ratio. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
October 28, 1981 
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TABLE I.- NOZZLE PRESSURE-ORIFICE LOCATIONS 0 


Configuration 1 

Configuration 2 

Configuration 3 

Configuration 4 

Configuration 5 

External 

Internal 

External 

Internal 

External 

Internal 

External 

J Internal 
! 

External 

Internal 

x/d m 

x/d m 

x/d m 

x/d m 

x/d m 

x/d m 

x/d m 

x/d Bl 

x/d m 

x / d m 



single external-orifice row is located along top center line of nozzle afterbody. Single 
internal-orifice row is located 45° clockwise from top of nozzle, looking upstream. 

^ Subtract these values from external and internal coordinates to locate orifice from start of 
boattail (station 137.16 cm). 
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Figure 1.- Photograph showing installation of nozzle configuration 5 in test section of 

the Langley 16-Foot Transonic Tunnel. 
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Figure 3.- Photograph showing the five nozzle configurations tested 
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Nozzle external pressure distribution 
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Nozzle internal pressure distribution 
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Figure 8.- External and internal surface-pressure distributions 
of nozzle configuration 4. Origin is at start of boattail; 
location of throat is denoted by vertical line in nozzle cross- 
section contour; design nozzle pressure ratio is 6.23. 
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Figure 8.- Continued. 

















Nozzle external pressure distribution 
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Figure 10.- Variation of jet-off, boattail pressure-drag coefficient with 

free-stream Mach number. 
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Figure 11 


Continued 
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Configuration 1. Design nozzle pressure ratio is 4.25. 



Configuration 2. Design nozzle pressure ratio is 21.23. 
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(a) Configurations 1 and 2. 

Figure 25.- Comparison of experimentally and analytically determined variation of 
static thrust coefficient with nozzle pressure ratio. 



Configuration 3. Design nozzle pressure ratio is 5.03. 
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Configuration 4. Design nozzle pressure ratio is 6.23. 
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(b) Configurations 3 and 4. 


Figure 25 Continued. 



Configuration 5. Design nozzle pressure ratio is 10.64. 



(c) Configuration 5. 
Figure 25.- Concluded. 
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